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Kaon decay studies seeking new-physics (NP) effects in leptonic {Kn) or semileptonic (K13) 
decays are discussed. A unitarity test of the first row of the CKM mixing matrix is obtained 
from the KLOE precision measurements of K13 widths for K , Kl, and (unique to KLOE) 
Ks, complemented with the absolute branching ratio for the A M 2 decay. KLOE results lead 
to constraints for NP models and can probe possible charged Higgs exchange contribution in 
SM extensions with two Higgs doublets. The main focus in the present document is set on a 
new measurement of Rk = T(K e 2) /T(K,j.2) with an accuracy at the % level, aiming at finding 
evidence of deviations from the SM prediction induced by lepton-flavor violation NP effects. 



1 Introduction 

New precise measurements of K —* lvi(pf) (K12) and K — > 7^1^(7) (^3) decays can possibly shed 
light on new physics (NP). The first indication of the need of improving the present knowledge 
in this field was given by the 2004 version of the PDG: a deviation from unitarity of the CKM 
matrix was observed in the first row, amounting to more than two standard deviations^, 

A = 1 - Vl d - Vl - Vl h = 0.0043(16)v^(n W (1) 

This called for new precise determinations of the V us parameter of the CKM matrix, traditionally 
extracted from K13 decays using the following expression: 

ri (^e3( 7 ),M3(7)) = 1^*1 1287T 3 ^ W ^| 2 ^3,Ai3 i 1 + 5 e3,^3)> ( 2 ) 

where i indexes K° — » ir~ and K + — > ir° transitions for which Cf = 1 and 1/2, respectively, Gf 
is the Fermi constant, M is the appropriate kaon mass, and /Sew is a universal short-distance 
electroweak correction^. The 5 l term accounts for long-distance radiative corrections depend- 
ing on the meson charges and lepton masses and, for K^, for isospin-breaking effects. These 
corrections are presently known at the few-per-mil level^. The /+"°(0) form factor parametrizes 
the vector-current transition K° — » tt~ at zero momentum transfer t, while the dependence of 



vector and scalar form factors on t enter into the determination of the integrals I e 3 u3 of the 
Dalitz-plot density over the physical region. 

After four years of analysis of KLOE data, we present the most comprehensive set of results 
from a single experiment, including BR's for K e % and Ku3 decays for Kl® and .fT 1 * 1 ^, and 
the BR for K$ — > Trev^^i* (unique to KLOE); form factor slopes from analysis of Kl £ ^^ and 
Klu3 > lifetime measurements for if^^l anc j ^±111 the j(0 mass L12] Using the K$ lifetime 
from PDGE1 as the only input other than KLOE measurements, we obtain five results for the 
product / + (0)|Vu S |Q3] as shown in table [TJ The average of these has been obtained taking all 
correlations into account and it is /+(0) x \V US \ = 0.2157(6). As a comparison, using data from 
KLOE, KTeV, NA48, and ISTRA+ experiments, the world average^ j s 0.2166(5). From the 
KLOE result and using / + (0) = 0.9644(49) from the UKQCD/RBC collaboration^, we obtain 

\V US \ = 0.2237(13). (3) 

Using the world average^l V u d = 0.97418(26) from + — > + nuclear j3 decays, CKM unitarity 
can be seen to be satisfied: A = 9(8) x 10~ 4 . 

KLOE has provided the most precise determination of the 

brUSI, which can be linked 

to the ratio V us /V U( i via the following relational 
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The theoretical inputs are the form- factor ratio fx/ fn and the radiative corrections described 
by the factor C. We use f K /U = 1.189(7) from lattice calculations^ an d c = 0.9930(35)121, 
thus obtaining 

\V US /V ud \ = 0.2326(15). (4) 

From the KLOE results of Eqs. [3]and|31 and from the world-average value of V u d, a combined 
fit to V us and V uc i has been done. The result is shown in left panel of Fig. [TJ the fit x 2 is 
2.34 for one degree of freedom (13% probability) and the results are: \V US \ = 0.2249(10) and 
\V u d\ = 0.97417(26), with a correlation of 3%. From these, not only can we now state that the 
CKM unitary holds to within 10" 3 , A = 0.0004 ± 0.0005^ ± 0.0004^, but we can obtain 
severe constraints for many NP models. 



1.1 Unitarity and coupling-universality tests 

In the SM, unitarity of the weak couplings and gauge universality dictate: 



Gl (|K 



,d\ 2 + \Vus\ 2 ) = Gl (V* b negligible , (5) 



where is the decay constant obtained from the measurement of the \x lifetime^!. The above 
measurement of V^ s from KLOE inputs provides relevant tests for possible breaking of the CKM 
unitarity (A ^ 0) and/or of the coupling universality {Gf ^ G>). This can happen in some NP 
scenarios, some example of which we discuss below. 

NP might lead to exotic and still unobserved [i decays contributing to the /i lifetime. The 
resulting total BR for fi exotic modes equals the unitarity violation A. Some of these modes, 
such as — * e + V e i' f j i , are at present constrained to be less than ~ 1%, so that information from 
unitarity improves on that from direct searches by more than a factor of 10 ^^ ! 

The existence of additional heavy Z bosons would influence unitarity at the loop level entering 
in muon and charged current semileptonic decays differently^^: A = — 0.01Aln[r^/(r^ — 1)], 



where rz = mz'/mw and A is a model-dependent constant of order 1. In the case of 50(10) 
grand unification, A ~ 1.9 and a unitarity test from KLOE results yields Mz' > 750 GeV at 
95% of CL. In non-universal gauge interaction models, a tree-level contribution from Z' bosons 
appears, so that the unitarity test is sensitive to even larger masses^S. 

In supersymmetric extensions of the SM (SUSY), loops affect muon and semileptonic decays 
differently. Unitarity can constrain SUSY up to mass scales of the order of 0.5 TeV, depending 
on the extent of cancellation between squark and slepton effects^! 

Measurements of Ki% widths can be linked to new physics effects, too. The ratio of K„2 to 
7Tu2 decay widths might accept NP contributions from charged Higgs exchange EZBEl j n super- 
symmetric extensions of the SM with two Higgs doublets. In this scenario, the ratio V us /V u d 
extracted from Ku2, 7t7t2 should differ from that extracted from K13 and superallowed Fermi 
transitions ("0+"): 



V us (K l2 )V ud (0 + ) 
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where tan/3 is the ratio o f up - and down-Higgs vacuum expectation values, Mjj is the charged 
Higgs mass, and €q ~ 0.01^1 The KLOE result of Eq. H]can be translated into an exclusion plot 
in the plane tan (3 vs Mh (see right panel of Fig. [1]), showing that this analysis is complementary 
to and competitive with that'SS using the average BR(I? — > tv) = 1.73(35) x 10 4 of Babar and 
Belle measurements^!. 
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Table 1: Five determinations 
of /+ x \V US \ using the K s 
lifetime (from PDG) as the 
only input other than KLOE 
measurements. 



Figure 1: Left: The 1-cr fit result to V ud and V us is shown 
by the solid line ellipse, in agreement with the unitarity 
bound shown by the dashed line. Right: Excluded regions 
from analysis of decays K — > [iv (filled area) and B — > tu 
(hatched area). 



1.2 Test of lepton-flavor violation 

A significant effort has been devoted along the years to isolate signals from lepton flavor violating 
(LFV) transitions, which are forbidden or ultra-rare in the Standard Model (SM). The sensitivity 
to decays such as \i — > &y, \i — > eee, Kl — > fie(+-K 0, s), and others roughly improved by two orders 
of magnitude for each decade^H. No signal has been observed, thus ruling out SM extensions 
with LFV amplitudes with mediator masses below ~ 100 TeV. 

These results allowed the focus to be put on the detection of NP-LFV effects in loop ampli- 
tudes, by studying specific processes suppressed in the SM. In this field, a strong interest for a 
new measurement of the ratio Rk = T{K — > ev)/T(K — > \xv) has recently arisen, triggered by 
the work of Ref.^Zl. The SM prediction of Rk benefits from cancellation of hadronic uncertain- 
ties to a large extent and therefore can be calculated with high precision. Including radiative 
corrections, the total uncertainty is less than 0.5 per mil^^ 



R K = (2.477 ±0.001) x 10~ 5 . 



(6) 



Since the electronic channel is helicity-suppressed by the V — A structure of the charged weak 
current, Rk can receive contributions from physics beyond the SM, for example from multi-Higgs 
effects inducing an effective pseudoscalar interaction. It has been shown in Refj^that deviations 
from the SM of up to few percent on Rk are quite possible in minimal supersymmetric extensions 
of the SM and in particular should be dominated by lepton-flavor violating contributions with 
tauonic neutrinos emitted in the electron channel: 



where Mjj is the charged-Higgs mass, A^ 1 is the effective e-r coupling constant depending on 
MSSM parameters, and tan (3 is the ratio of the two vacuum expectation values. Note that the 
pseudoscalar constant fx cancels in R^. 

In order to compare with the SM prediction at this level of accuracy, one has to treat carefully 
the effect of radiative corrections, which contribute to nearly half the width. In particular, 
the SM prediction of Eq. [7] is made considering all photons emitted by the process of internal 
bremsstrahlung (IB) while ignoring any contribution from structure-dependent direct emission 
(DE). Of course both processes contribute, so in the analysis DE is considered as a background 
which can be distinguished from the IB width by means of a different photon energy spectrum. 

Two experiments are participating in the challenge to push the error on Rk from the present 
6% down to less than 1%. In 2007, KLOE and NA48/2 announced preliminary results'^ with 
errors ranging from 2% to 3%. Moreover, the new NA62 collaboration collected more than 
100 000 K e 2 events in a dedicated run of the NA48 detector, aiming at reaching an accuracy of 
few per mil on Rk^^- 

2 Measuring R K at KLOE 

DA<I>NE, the Frascati factory, is an e + e~ collider working at y/s ~ ms ~ 1.02 GeV. <j> mesons 
are produced, essentially at rest, with a visible cross section of ~ 3.1 fib and decay into K + K~ 
pairs with a BR of ~ 49%. 

Kaons get a momentum of ~ 100 MeV/c which translates into a low speed, 0k ~ 0.2. K + 
and K~ decay with a mean length of X± ~ 90 cm and can be distinguished from their decays 
in flight to one of the two-body final states [iv or tttt . 

The kaon pairs from (j) decay are produced in a pure J PC = 1 quantum state, so that 
observation of a K + in an event signals, or tags, the presence of a K~ and vice versa; highly 
pure and nearly monochromatic beams can thus be obtained and exploited to achieve high 
precision in the measurement of absolute BR's. 

The analysis of kaon decays is performed with the KLOE detector, consisting essentially of 
a drift chamber, DCH, surrounded by an electromagnetic calorimeter, EMC. A superconducting 
coil provides a 0.52 T magnetic field. The DCIP^is a cylinder of 4 m in diameter and 3.3 m in 
length, which constitutes a fiducial volume for decays extending for ~ 1A±. The momentum 
resolution for tracks at large polar angle is cr p /p < 0.4%. The cm. momenta reconstructed from 
identification of 1-prong — > fiv, tttt decay vertices in the DC peak around the expected 
values with a resolution of 1-1.5 MeV, thus allowing clean and efficient tagging. 

The EMC is a lead/scintillating-fiber sampling calorimeter^] consisting of a barrel and two 
endcaps, with good energy resolution, cte/E ~ 5.7%/v E(GeV), and excellent time resolution, 
a T = 54 ps/VE(GeV) 50 ps. 

In early 2006, the KLOE experiment completed data taking, having collected ~ 2.5 fb _1 
of integrated luminosity at the 4> peak, corresponding to ~3.6 billion K + K~ pairs. Using the 
present KLOE dataset, a measurement of Rk with an accuracy of about 1 % has been performed. 





Given the decay length of ~90 cm, the selection of one-prong K^ 1 decays in the DC 
required to tag has an efficiency smaller than 50%. In order to keep the statistical uncertainty 
on the number of K — > ev counts below 1%, a "direct search" for K —* ev and K —> fj,v decays 
is performed, without tagging. Since the wanted observable is a ratio of BR's for two channels 
with similar topology and kinematics, one expects to benefit from some cancellation of the 
uncertainties on tracking, vertexing, and kinematic identification efficiencies. Small deviations 
in the efficiency due to the different masses of e's and //'s will be evaluated using MC. 

Selection starts by requiring a kaon track decaying in a DC fiducial volume (FV) with labo- 
ratory momentum between 70 and 130 MeV, and a secondary track of relatively high momentum 
(above 180 MeV). The FV is defined as a cylinder parallel to the beam axis with length of 80 cm, 
and inner and outer radii of 40 and 150 cm, respectively. Quality cuts are applied to ensure 
good track fits. 

A powerful kinematic variable used to distinguish K — > ev and K — > \iv decays from the 
background is calculated from the track momenta of the kaon and the secondary particle: as- 
suming My = 0, the squared mass of the secondary particle (M± ) is evaluated. The distribution 
of M[g p is shown in Fig. [2] for MC events before and after quality cuts are applied. The selection 
applied is enough for clean identification of a K — > [iv sample, while further rejection is needed 
in order to identify K — > ev events: the background, which is dominated by badly reconstructed 
K — ► p,v events, is ~10 times more frequent than the signal in the region around M c 2 . 
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Figure 2: MC distribution of M? 
before (left) and after (right) quality 
cuts are applied. Shaded histogram: 
K — > ev events. Open histograms: 
background. In MC, Rk is set to the 
SM value. 
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Figure 3: Left: cell distribution for 200 MeV e 
(top) and fjL (bottom) from two selected events from 
K L -> ir£v. Right: Distribution of NN output, NN, 
for electrons of a Kl — > nev sample from data (black 
histogram) and MC (red histogram). 



Information from the EMC is used to improve background rejection. To this purpose, we 
extrapolate the secondary track to the EMC surface and associate it to a nearby EMC cluster. 
For electrons, the associated cluster is close to the EMC surface and the cluster energy E c \ is 
a measurement of the particle momentum p ey ±, so that E c \/p eX t peaks around 1. For muons, 
clusters tend to be more in depth in the EMC and E c \/p ext tends to be smaller than 1, since 
only the kinetic energy is visible in the EMC. Electron clusters can also be distinguished from \i 
(or 7r) clusters, since electrons shower and deposit their energy mainly in the first plane of EMC, 
while muons behave like minimum ionizing particles in the first plane and deposit a sizable 
fraction of their kinetic energy from the third plane onward, when they are slowed down to rest 
(Bragg's peak), see left panel of Fig. [3l Particle identification has been therefore based on the 
asymmetry of energy deposits between the first and the next-to-first planes, on the spread of 
energy deposits on each plane, on the position of the plane with the maximum energy, and on the 
asymmetry of energy deposits between the last and the next-to-last planes. All information are 
combined with neural network (NN) trained on Kl — > irttv data, taking into account variations 
of the EMC response with momentum and impact angle on the calorimeter. The distribution of 
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Figure 4: Fit projections onto the My axis for two 
slices in NN output, NN > 0.98 and NN < 0.98, giv- 
ing enhanced values of signal and background contri- 
butions, respectively. 




Figure 5: 95%-CL excluded regions in the 
plane tan /3-charged Higgs mass for A 3 ^ 1 = 
10" 4 , 5 x 10~ 3 , 10~ 3 . 



the NN output, NN, for an independent Kl — ► neu sample is shown in the right panel of Fig. [3] 
for data and Monte Carlo (MC). Additional separation has been obtained using time of flight 
information. 

The number of K — > ev{pf) is determined with a binned likelihood fit to the two-dimensional 
NN vs M[g p distribution. Distribution shapes for signal and background, other sources being 
negligible, are taken from MC; the normalization factors for the two components are the only fit 
parameters. In the fit region, a small fraction of K — > ev{^) events is due to the direct-emission 
structure-dependent component (DE): the value of this contamination, /go, is fixed in the fit to 
the expectation from simulation. This assumption has been evaluated by performing a dedicated 
measurement of SD, which yielded as a by-product a determination of fsD with a 4% accuracy. 
This implies a systematic error on K e 2 counts of 0.2%, as obtained by repeating the fit with 
values of /sd varied within its uncertainty. 

In the fit region, we count 7064±102 K + — > e + v{~f) and 6750±101 K~ — > e~&('y) events. 
Fig. [5] shows the sum of fit results for K + and K~ projected onto the My axis in a signal- 
(NN > 0.98) and a background- (NN < 0.98) enhanced region. 

To assess the uncertainty on the Rk measurement arising from limited knowledge of the 
momentum resolution we have examined the agreement between the Mf ep distributions for data 
and MC in the K^2 region. For the NN distribution, the EMC response at the cell level has 
been tuned by comparing data and MC samples. In order to evaluate the systematic error 
associated with these procedures, we studied the result variation with different fit range values, 
corresponding to a change for the overall K e 2 purity from ~ 75% to ~ 10%. The results are 
stable within statistical fluctuations. A systematic uncertainty of 0.3% for Rk is derived "a la 
PDC'^^by scaling the uncorrelated errors so that the reduced x 2 value of results is 1. 

The number of events in the same data set is extracted from a fit to the My distribution. 
The fraction of background events under the muon peak is estimated from MC to be < 0.1%. 
We count 2.878 x 10 8 (2.742 x 10 s ) K+ 2 (K~ 2 ) events. Difference in K+ and K~ counting is 
ascribed to K~ nuclear interactions in the material traversed. 

The ratio of K e 2 to efficiency is evaluated with MC and corrected for data-to-MC ratios 
using control samples. To check the corrections applied we also measured R3 = BR(K e ^) /BR^i^), 
in the same data sample and by using the same methods for the evaluation of the efficiency as 
for the Rk analysis. We found R3 = 1.507(5) and R3 = 1.510(6), for K + and K~ respectively. 
These are in agreement within a remarkable accuracy with the expectation^^ from world-average 
form-factor slope measurements, R3 = 1.506(3). 



3 Rk result and interpretation 



The final result is R K = (2.493 ± 0.025 ± 0.019) x 10~ 5 . The 1.1% fractional statistical error 
has contributions from signal count fluctuation (0.85%) and background subtraction. The 0.8% 
systematic error has a relevant contribution (0.6%) from the statistics of the control samples 
used to evaluate corrections to the MC. The result does not depend on K charge: quoting olny 
the uncorrelated errors, R K {K + ) = 2.496(37)10~ 5 and R K {K~) = 2.490(38)10" 5 . 

The result in agreement with SM prediction of Eq. [6l Including the new KLOE result, the 
world average reaches an accuracy at the % level: Rk = 2.468(25) x 10~ 5 . In the framework of 
MSSM with LFV couplings, the Rk value can be used to set constraints in the space of relevant 
paremeters (see eq. [7]). The regions excluded at 95% C.L. in the plane tan /3-charged Higgs mass 
are shown in Fig. [5] for different values of the effective LFV coupling . 
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